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Abstract 

The paper revisits the compressible Kraichnan model of turbulent advection in order to derive 
explicit quantitative relations between scaling exponents and Lagrangian particle configuration 
geometry. 



1. Introduction 

Physical and numerical experiments have typically access to, at least, two kinds of observ- 
ables of the state of a turbulent Newtonian fluid. These observables are, on the one hand, Eulerian 
or Lagrangian statistical indicators such as universal scaling properties of correlation and struc- 
ture functions. On the other hand, geometrical indicators like shapes formed by configurations 
privileged by the dynamics of Lagrangian or inertial particles can be observed. A series of re- 
markable results (see e.g. [1,2] and references therein) have in recent years started to unveil 
the nature of the link between statistical indicators and geometrical properties. In particular, [3] 
gave clear evidence that tracer particles advected by the 2d Navier-Stokes velocity field spend 
anomalously long times in degenerate geometries characterized by strong clustering and that 
the phenomenon has a quantitative counterpart in the existence of special functions of particle 
configurations which are on average preserved by flow. Furthermore, it was shown in [4] for den- 
sities transported by compressible turbulent velocity fields that statistical conservation laws de- 
termine, at least for integer values of the mass, the multifractal spectrum of the attractor towards 
which Lagrangian trajectories converge. On the background of these results was the discovery 
of anomalous scaling [5, 6, 7] in a stylized model of passive advection, the Kraichnan model 
[8]. There it was possible to show in a mathematically controlled fashion [9], that the universal 
statistical properties of the model in the inertial and, under additional hypotheses [10], decay 
ranges are determined by statistical conservation laws of the Lagrangian dynamics (referred to 
as zero modes). Zero modes have subsequently become the object of extensive investigations 
and analytical, mostly perturbative, and numerical expressions of their scaling exponents have 
been derived for the Kraichnan model and its generalizations (see e.g. [11, 12] and references 
therein). To the best of my knowledge, however, quantitative expressions of the scaling exponents 
explicitly relating them to geometrical properties of Lagrangian particle configurations appeared 
only recently in [13]. The scope of the present contribution is to illustrate the new method of 
calculation introduced in [13] in the slightly more general case of the compressible Kraichnan 
model [14]. The main result of the paper, derived in sections 4 and 5, is the classification of zero 
modes, irreducible and reducible, in terms of Gel'fand-Zetlin patterns (see e.g. [15]) associated 
to quadratic Casimir operators of classical groups. It must be clearly stated that this classification 
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is peiturbative but it has nevertheless the merit to establish an explicit link to Lagrangian particle 
geometry which is discussed in more details for four point correlations in section 6. Sections 2 
and 3 recall respectively the defining properties of the Kraichnan model and basic facts about the 
structure of the solution of the Hopf equations. In the conclusions I discuss the use of the new 
method and make some speculations on its significance for the Navier-Stokes equation as well 
as for other statistical field theories. 



2. The compressible Kraichnan model 

Passive turbulent advection by a compressible velocity field in M*^ encompasses (see e.g. 
[14]) the density evolution equation 

dtP + da,-{vp)='^dlp + f (1) 

(describing e.g. the distribution of particles floating on the surface of the fluid in the limit of 
vanishing Stokes number) and the tracer evolution equation 

dtO + v da^e = '^d'^e + g (2) 

(describing e.g. the fluid temperature). At finite molecular viscosity k, the Lagrangian dynamics 
underlying (1) and (2) are distinct, respectively forward and backward in time. In the inviscid 
limit and in the case of advection by a turbulent field stylized by a stationary, time-decorrelated 
Gaussian ensemble they become formally adjoint [14] and can be discussed in parallel. These 
are the working hypotheses of the paper. In particular, the velocity field in (1), (2) is modeled by 
the realizations of Kraichnan's zero average Gaussian compressible ensemble: 

~< v''{xuti)v^{x2M) y= ^(ti2)£>"^(xi2,m) (3a) 



D<^^ix;m) = DoH2-0 / (0^/^ ' P)'"" ' ' ' ^ pf-^} X [f) (3b) 

where Xi2 '■— Xi — X2, ti2 := <i — h, ^ & [0, 2] is the roughness degree of the velocity field, 
p G [0, 1] is the degree of compressibility and x some non-universal infra-red cut-off function 
normalized to the unity for vanishing inverse integral scale m. The force fields /, g (Gaussian, 
zero-average and decorrelated in time) compensate the dissipation, in the inviscid limit only due 
to the eddy diffusivity generated by the velocity field, in order to let the system attain a Galilean 
invariant steady state. Consistence with hydrodynamics imposes to interpret (1), (2) advectedby 
the Kraichnan ensemble as stochastic partial differential equations in the sense of Stratonovich. 
The role of the eddy diffusivity appears upon converting (1), (2) into Ito form [11,12]. This latter 
is for the density 

dtP + d^-{vp) = ^:^^d'p + f (4) 

and for the tracer 

dte + vd^e = ^:^^d''e + g (5) 
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The eddy diffusivity in the above equations is specified by 



which means that Taylor's fomiula [16] becomes exact in the Kraichnan model. 

3. Hopf equations and martingales 

A straightforward application of Ito calculus to (2) shows that the equal time correlation C„ 
of n tracer fields in R'' satisfies the Hopf equation (see e.g. [9, 11]) 

(dt + M^^) C„ {X; t) = J'n {X; t) (7) 

n 

= - 2 E ^ (Xi,- , m)d,. d^, (8) 

ij 

and J>j an effective forcing term fully specified by the correlation of g and by correlation func- 
tions Cn' with n' < n. The operator Mn can be also regarded as the generator of a multiphcative 
diffusion process the transition probabiUty Vn whereof obeys 

{dt+M'^^^)Vn{X,X'-t-t')=0 (9a) 



with 



Vim.Vn{X,X';t-t') = d^'^^'HX - X') (9b) 
tit' ^ J ^ I 

The Lagrangian interpretation of Vn is that it is the probability density to find a cluster of 
n inertial particles in X' = (aj^ , . . . , a;^) at time t' conditioned upon the event that they reach 
X = {x\, . . . , Xn) at a later time t. Analogously, the equal time n-point correlation of the 
density field satisfies an Hopf equation formally adjoint to (7). In such a case Vn can be regarded 
as the solution of the Fokker-Planck equation 

(dt + M^;p,^)Vn{X,X';t-t') =0 (10) 

with 

1 " 

M^x'^ = -,^Y.^.t^^'D'^^{xi,,m) (11) 

ij 

Vn acquires the interpretation of the probability density for a cluster of n inertial particle to arrive 
in X' = (x[, . . . , x'^) at time t conditioned upon the event that they are in X = (xi, . . . , x„) 
at time t' < t. The different interpretations of (9a) and (10) stem from the distinct Lagrangian 
dynamics underlying (1), (2) at finite molecular viscosity [14]. Whenever the kernel Vn acts 
on translational invariant functions, the arithmetic average of the Lagrangian particle positions is 
integrated out. The reduction defines a d„ = {n—1) d-dimensional subspace of M" the ttansla- 
tional invariant sector of the theory. The projection P„ of Vn on the subspace has a well defined 
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scale invariant limit for m tending to zero [9]. The scaling dimensions of P„ corresponding to a 
linear rescaling of spatial variables (i.e. in units d^:) are 

dt = (2-e)d^ & dp„ = -(n- l)dd^ = -d„d^ (12) 

Qualitative analysis [9] yields for P„ the asymptotic expansion: 

i y2Tn'Pi-iiY)ipi-i(Y',t-t') fory<F' 
Pn{Y,Y',t-t') = \ ^^="^'''^ ^^''^^ ' - , Y,Y'eR''"(13) 

Throughout the manuscript, Ha denotes the Heaviside function with normalization Ha{0) = a. 
The meaning of the expansion is well illustrated by the £, = limit. In such a case 

ywjj- —Jvi^ — 2 ~ 2 

and both Vn and the reduced propagator P„ are Gaussian. Self-adjointness of the Laplacian also 
implies 

fpi = = V'i;2 & (t>i = = ^i;2 (14) 

It is expedient to couch the sum over i into a sum over {j,l,k),j being the degree of homogeneity 
of an harmonic polynomial Hji, I the numbers specifying an S'0(d)-adapted representation 
of hyperspherical harmonics of SO{dn) [17] and 

MY) ■■= ^-^^^^ V - y^^-^y,iY/y) 

22fcr(fc + i)r( '^"+^/+^" ) 22fcr(fc + i)r( ^"+^/+^" ) 

for y = ll'Fll. The so defined 0jj/c;i's satisfy the "tower" relations 

d'^(l)jik{Y) = cl)jik-i{Y), cl>ji_^{Y) = (16) 

Under the same conventions, the identity 

2 I d„,+2 j-2 / „2 \ 

MY, t) = 2n^ ^^^^ 

holds true for the generalized Laguerre polynomial of degree k. An alternative useful rep- 
resentation is 

(vn-.^^ {Y,Kt\\jlk)e-^^ 
%lk(Y,t) = 2 7r 2 — — — (18) 

with {Y,>ct\ \jlk) the wave function of the {j,l,k) eigenstate of the isotropic quantum har- 
monic oscillator of unit mass and frequency in IR''"[18]. Finally, the ipjik{Y, t)'s satisfy the 
"tower" relations 

dti^MY,t) = ^d^i^jik{Y,t) = -ijjik+i(Y,t) (19) 
Some remarks are in order. 



i The identities (16), (19) are particular case of consistence conditions of the expansion (13) 
holding for generic ^ |9, 11]. 

a By (18) and the orthonormality of quantum eigenstates it readily follows that the harmonic 
polynomial is L^(]R'*'')-orthogonal to the V'jJfe's for any fc > 1. The normalization 
of the (15)'s is chosen such that they are exactly replicated whenever averaged over the 
transition probability P„. In other words, they are martingales of the Wiener process [ 19]. 

Hi The time integral 

[dti,%[Y,t)^y,^l^^ (20) 

defines for k equal zero only local martingales (i.e. harmonic functions). The singularity 
for y J, 0, obstructs (20) absolute integrability and therefore the martingale property. This 
fact has important consequences for perturbative investigation of anomalous scaling in the 
decay range (scales larger than the integral scale of the forcing for infinite integral scale of 
the velocity field [13]). 

4. Scaling, martingales of the diffusion process and perturbation theory 

The term martingales [19] denotes functions of a Markov process the shape of which, on 
average, is preserved by the time evolution: 

j (fxT{x, s) K [x, s\x', s') = T{x\ s') (21) 

The kernel if in (21) is a probability density for x at time s conditioned upon the occurrence of 
x' at time s' . For forward (backward) dynamics s > s' (s' > s). Clearly, a time independent 
T{x) satisfying (21) is on average a conserved quantity of the diffusion process described by the 
transition probability K. Furthermore, (21) requires only integrability of J^{x) with respect to 
the transition probability. The functional space where statistical conservation laws can be sought 
is therefore larger than that of the physical solutions of the Hopf equations. Turning the attention 
to the dual action of K, the role of constant martingales is played by stationary pseudo-measure 
Q satisfying 

j (fx' K {x, s\x', s') g{x') = g{x) (22) 

The relevance of these considerations for the Kraichnan model [9, 11] is that the expansion 
(13) generically comprises constant martingales Zi-iand stationary pseudo-measures Zi.2- They 
satisfy the zero mode equations 

7W(")2:,;i = & A^'"^^Z,;2 = (23) 

and asymptotically dominate in the inertial range the scaling of respectively the tracer and density 
correlation functions. Observation Hi at the end of the previous section evinces that the concept 
of statistical conservation law (i.e. martingale or stationary pseudo-measure) is stronger than 
the one of zero mode as the (23)'s provide only necessary conditions. For small ^, small scale 
zero modes are deformations of harmonic polynomials which are statistical conservation laws of 
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the Wiener process. In such a case the identification of zero modes with statistical conservation 
laws holds true. The fact suggests to use (21), (22) to device a perturbative scheme to determine 
the scaling dimensions Q-r of zero modes Zi-r- Instead of solving perturbatively (23) as in 
[5, 20], it is possible to read the dependence of scaling dimensions upon ^ by looking at how 
the transition probability density deforms statistical conservation laws of the Gaussian theory. 
If statistical conservation laws exist for small but finite values of the roughness parameter, the 
deformation introduces a logarithmic time dependent prefactor in front of harmonic polynomials. 
Conservation laws at finite ^ must cancel this time dependence by containing order by order in ^ 
suitable logarithmic counter-terms in the spatial variables. The constant prefactor of these terms 
yields the coefficients of the Taylor expansion of the (i;r's,- In formulas, the claim is that choosing 
the state numbers J = (j, I) so that the action of the transition probability is within logarithmic 
accuracy diagonal then 

Zj = e"*-^" o f — I Hj + subleading or amplitude terms (24) 

Since the perturbative expansion of a statistical conservation law Zi cannot contain any time t 
dependence, the scaUng dimensions of Z and Z must satisfy 

d2 = dj + [Cj(0-O(0)] (25) 



5. Calculation of scaling dimensions 

As solution of (9a) at leading order in perturbation theory in ^ the transition probabiUty 
density Vn is 

e-^" - ^^^A^e^^" + f + Oie) (26) 

The subscript n/{ij) betokens dependence upon all n particles coordinates but the i-th and the 
j-th particle. 



g-tM2 (27) 

€=0 



is the first order correction to the two-points response function. 

5.7. Diagrammatic expression o/r'^^ 

R2'^ can be computed by standard diagrammatic techniques of statistical field theory (see 
e.g. [12]). Denoting for t > s 

^ ^ 
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and derivatives by a line perpendicular to response lines, r'^^ is given by 



R^2\xi,Xj,t\yi,y^,s) 



• AA/H^AA* 

, ^, • — v\aJh^w» , , 



with 



PlaP20 (^l-e-*^[Po+Po-(Pi-P2 



and 

ro := Xi - Xj , n := Xi-yi, r2 := - y^- (30) 
A (gradient) expansion in the momenta p^, pj recasts the integral (28) into the form 



R^2\xi,Xj,t- s\yi,yj,0) = '^vl^^l^{xi,Xj,t- s)e * \xi,Xj\yi,yj) (31) 

fe=o 



(k) 

The vertices V^^^ ^ 's are homogeneous normal ordered differential operators of degree k. Scaling 

analysis requires to focus on the terms [V^^f^''^] in the V^^'^^'?, proportional to \n{xt)^^ . Only the 
explicit expression of the first few of them is needed. The reason is the following. Upon inserting 
(31) into (26), the transition probabiUty reduces to 



= e^A„_^^^^^_^^O^lM[e^-"+0(a (32) 

i^j k 

Harmonic polynomials Hji behave as eigenvectors of unit eigenvalue of the propagator: 

i^j k 

Suppose that T-L^i is particle-permutation and translational invariant, as it is expected for inertial 
range tracer correlations. Suppose furthermore that 'Hji contains powers of (components of) the 
variable Xi not lower than k' < j i.e that it has minimal homogeneity degree k'. Then such 
polynomial can only be reproduced by operators [V^^j'^^' ^ ] with k < k'. The conclusion is that 
scaling analysis at order ^ requires in such a case the spectrum of the homogeneous differential 
operator 

<.^) = ^EE<.)](-^.-i) (33) 

k=0 i^j 

Physically, special interest have conserved quantities dominating the inertial range scaling of 
structure functions of the tracer. Symmetry and analyticity then restrict the corresponding scaling 
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(k) 

analysis to the diagonalization of all the 0^^\'s for fc < 2. At this stage it is worth observing 
that the argument above also applies to the scaling analysis of the density field. The action of the 
transition probability on conservation laws of the dual space is 

e-t^l^^ = nj + ^J2Il ^i^t) + 0{e) (34) 

fc 

where the operators V^^^^ ^'s are obtained by replacing in (29) r!^^ with 



^^'^^Po.Pi.Pa,*) := , r \^ -T^ -DJ^^iPo) (35) 



5.2. Scaling analysis for the tracer 

The only vertex contributing to (33) is 



A ...... 

{2'n:Y xq^ 



V^l^(xux,,t) = / ^ e---^- A ^ "^'^ (36) 



By (24) in order to determine scaUng it is sufficient to extricate the logarithmic asymptotics of 
(36) for 



|2 



Kt 



<. 1 (37) 



The technical details of the calculation of the integral have no conceptual relevance and are out- 
lined in Appendix A. After some algebra (Appendix C), it tums out that translational invariant 
harmonic polynomials of minimal homogeneity degree two can specify the limit of vanishing 
^ of zero modes of (9a) if they coincide with the harmonic component of the homogeneous 
polynomials diagonaUzing the operator 

.42) _^ [V(^^o)](a;^,a;j) _ rf+l-2p f^(2,„) + (rf - 2) p ^(2,„) 
^{tr.)-2^ 2 ~ 2(d + 2)(d-l) r^O(<') d+l-2p ^f^("-i) 



^_{[d+{d-2)p\{d+\-n)[En+d{n- 1)] 
2d \ 



2d \ {d + 2){d-l){n-l) 



-d)^ Er, 



+ piEn -d)\ E„ (38) 



where C^^^^^ and are respectively the representations of the quadratic Casimir of 

SO{d) and SU (n — 1) on the space of functions of n-particles in c?-spatial dimensions (see Ap- 
pendix B). As c'^Q(liy ^sj/"!-!) commute, the spectrum of oj^^ ^ diagonalized on the space of 
translational and permutation invariant homogeneous polynomials of degree j can be classified in 
terms of the Gelfand-Zetlin pattems (see e.g. [15]). In essence, denoting by a = [ai, . . . , a„_2] 
n — 2 non-negative integers satisfying 

n-2 

ai = J & Oj > ai' for i> i' (39) 



the eigenvalues of C^^(l^_i) read 



n-2 



Xsuin-1) (a) = ^ ai (a, - 2 i) + (40) 



whilst the Casimir of C^^"^^ is 



Aso(d)W=^(^ + rf-2) (41) 

for i the total (hyper)-angular momentum. The conclusion is that permutation and translational 
invariant zero modes dominating the inertial range asymptotics of the n-point tracer correlation 
functions scale with exponents 



Ci^{([ai,...,a„_2],£)=n 

ad+l-2p) [ d+(d-2)p 
+ 2(d + 2)(d-l) r^^ + "^ d+l-2p 



"n-2 
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.2 



ai{ai+n-2i) - 

n — 1 



-2d\ (d + 2)(d-l)(n-l) + ) (42) 

Of particular interest is the irreducible SO{d) isotropic zero mode governing the scaling of the 
structure function. This latter corresponds to the pattern of highest symmetry [n, 0, . . . , 0] and 
scales with exponent 

CS ([n, 0, . . . , 0], 0) = n - ^ M^ + n + 2^(n-2)p] ^ ^^^.^ ^^3^ 

a result first obtained in [21] using operator product expansion (see also [14] for discussion). 

5.3. Scaling analysis for the density 

The adjoint action of the transition probability density brings about four interaction vertices. 
The first does not bring about any differential operation on harmonic polynomials: 

,„, f ra MS (l-e-""') D'f{g;m) 

v;2.,(««.*):=-/(|;^«--"— ^ ;^ («) 

The other two, after use of the identities 

D''^{q,m)qc,qp = -pD<'^{q,m) (45) 

and 

D^f'{q,m)q^pp = -pD^^{q,m)q-p (46) 

reduce to 



and 



Vf2(x„x„t) = ^J ^e"'-^^e-*«'' D« (g,m) [q . {d^^ - d^,)]' 
The fourth vertex is V^^f^ = vj^f^^^ . As (47) and (48) are not logarithmic in 



[Vl2.)]{xi,x^) = [vll']]{xuxj)=0 



the only new contribution relevant for scaling analysis comes from V^'^l ^ and is 



(de.) 



[vf2,](-.-) = -^ 



(48) 



(49) 



(50) 



(36) indicates the existence of zero modes based at harmonic polynomials of zero degree in the 
limit C 4 0. For these zero modes the operator 



O 



(0) 
(de.) 



is diagonal and predicts the scaling dimension [21, 14] 

.pdn{n — 1) 



Cn;: 



+o{e) 



(51) 



(52) 



The necessary condition that harmonic polynomials of minimal degree two need to satisfy in 
order to specify the limit of vanishing ^ of zero modes is to specify the harmonic component of 
the homogeneous polynomials diagonalizing the homogeneous operator 



0(2) 



>dn (n — 1) 



+ 



,(2,n) + (d - 2) p^(2,n) 



2 ■ 2(d + 2)(rf-l) d+l-2p ^^("-i) 

1 nd+{d-2)p]{d+l-n)[Er, + d{n-l)] 



2d 



p{En -d) \ (53) 



(d + 2) 

Again, Gel'fand-Zetlin patterns satisfying (39) provide a natural way to classify the spectrum of 

\de.y 



0)J ^ . The corresponding spectrum of scahng dimensions is 



pdn{n-l) {d+l-2p)l{i + d-2) 



2{d + 2) 



d+{d-2)p 
2{d + 2) (d-l) 



'n-2 



^ ai{ai +n-2i) 



J-1 



6. Discussion 

The general expressions (42), (54) are fairly comphcated. It is convenient to discuss their 
significance in the particular non-trivial case n = 4. This is the case of zero modes contributing 
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to fourth order structure and correlation functions. Restricting the focus on isotropic zero modes, 

l(2) Ai(2) 
'itr.y '^{de.) 



(2) (2) 

the Gelf'and-Zethn patterns diagonahzing O)'^, O),' ^ select two harmonic polynomials of 



degree four 

-2 ^ {xi - Xj)'^{xi - Xj')'^ +6 ^ {xi- Xi')'^{xj - Xj')'^ (55) 

{{i,j},{hj'}} {{i,i'},{j,f}} 

and 

^([2,2],0)cx-(2(i+l) ^(xi-XiO^ + (rf + 2) Yl {xi-Xjf{xi-x,,f (56) 

{{i,3},{h3'}} 

where the pairs {i, i'} and { j, /} are assumed different, as well as the pairs {i, j'} and { j, i'}. 
(55), (56) are the generalization to d-dimensions of the expressions given in [5] for d = 3. They 

(2) (2) 

specify the harmonic component of the homogeneous polynomials diagonalizing O^^^ y O^^^ y 
They differ from the full eigenstate by linear combinations of slow modes i.e. polynomials of 
the form <j>f,i,j^f's with / < j which are L^-orthogonal to the harmonic part on the hyper- 
sphere (S''"^^ in the general case). In this sense, harmonic polynomials are rephcated by the 

A second relevant aspect [10, 22, 1 3] on which to pay attention is represented by the relations 
interweaving small and large scale scaling exponents. Consider the Green function associated to 
the translation invariant transition probability density P„. The asymptotic expansion 

^n'^^^')=t\t^^S^ IIy'^ MY):=fdt^,Ay,t) (57) 
dictates 

+ = d^^,, + d0„,, = [2 - e - d(n - 1)] d^ (58) 

Bearing in mind the physical interpretation of the adjoint action of P„ these relations state that 
the knowledge of the scaling dimensions of inertial range zero modes of the tracer field entails 
that of large scale zero modes of the density field and vice- versa. For example, denoting with 
superscripts (s) and (/) respectively small and large scale exponents, the isotropic irreducible 
and reducible zero modes of the tracer four point correlation function give 

Clf)([4,0],0) = 4-2^i^±g^ + O(e2) 

Cg([2,2],0) = 4-2C^+O(a 

^ d!k[4,0],0) = -(3d + 2) + e^^±f^ + O(e2) 

Cl|k[2, 2], 0) = -(3 d + 2) + ^^^ip + 0{e) 
whilst the corresponding object for the density yield 

Ci:j([4,0],0) =4-2g '^+^+^''f+f +0(^^) 
'(*)cr9 91 - /I _ 9 , ^/r2\ 



(59) 



2], 0) = 4 - 2 ^ '^-^-^'.^^'^"-^n + 

Cl|l([4, 0], 0) = -(3 d + 2) + g W2p[4+3W)] ^ ^(^2) 

^ Ci!l([2,2],0) = -(3d + 2)+e^^^2±^4t±^^^^^ + O(^2) 



(60) 
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As expected, in the incompressible limit of vanishing p the relations become self-dual. Large 
scale zero modes may become manifest in the power law decay of certain statistical indicators 
at large point separations if the integral scale of the velocity field is much larger than the one of 
the forcing [22, 13]. A final observation is that a direct determination of scaling dimensions of 
large scale zero modes is naturally achieved by studying the time dependent martingales tpj i o 
rather their stationary counterparts i o which are not martingales. Within the approximation 
(37), straightforward algebra yields 

j d'^-Y'Pn{Y,Y',t-s)i)nlfi{Y',s) -=%nlAY,t) + 

^ t)[-l In(v^) + + In -^^O^^^H^i | +••• + 0(^^X61) 

where t = t — s and dots stand for slow modes. Scaling prediction for time-dependent martin- 
gales can be then read from the prefactor of the logarithmic counter-terms needed to compensate 
the the ln(x-T)~^ dependence in (61). 

7. Conclusions 

The technique illustrated is strongly reminiscent of the operator product expansion [23] pre- 
viously also used to perform systematic calculations in the Kraichnan model (see e.g. [24, 21, 
25, 12] and references therein). The differences are that the role of momentum cut-off's is played 
by the time and that instead of bases of abstract operator valued fields (composite operators) it 
makes use of a basis of harmonic polynomials or shapes complete on the S''""^ hypersphere. 
These objects have a direct geometrical meaning and can be used in numerical simulations close 
to controlled or phenomenological Gaussian limits to test scaling properties as it was done in 
[13]. A final observation is that, at least in principle, the technique does not require linearity 
of the statistical field theory in order to be appUcable. The reason is the following. In the ther- 
modynamic formalism of field theory [23], connected n-point correlation functions PF^") are 
reconstructed by solving Hopf-like equations of the form 

t/(2) vk(") = J-„(Ty("-i\ VK^i)) (62) 

where U^"^^ stands for the set of proper vertices of order two acting on n-point connected cor- 
relation functions. The functionals J'„'s in (62) are specified by functional derivatives of the 
Legendre transform connecting the free energy W (i.e. the generating function of connected cor- 
relations) to the thermodynamic potential (i.e. the generating function of proper vertices). Thus 
the J>,'s depend upon aU proper vertices of order less or equal to n. Once these latter ones are 
given, (62) defines a solvable hierarchy for the connected correlations. Thinking of the ^/'■"^'s as 
pseudo-differential operators, the generalization of the method presented in this paper consists in 
probing the scaling of "martingales" associated to the generalized propagators U^'^'>~^ acting on 
n-particle functions. It must be, however, emphasized that deriving expressions of proper ver- 
tices is in general a very non-trivial task that can be most often carried out only when perturbative 
methods are applicable. 
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Appendix A. Outline of the evaluation of the integrals 



A convenient method to remove cut-offs from integrals associated to Feynman-diagram is to 
take the MelUn transform with respect to the cut-off (see e.g. [23]). For the velocity correlation 
this means 



{^,rn)- J ^^^^^e ^^^^ (A.l) 



with n(p, p) = {l- p)S"f^ - (1 -^dp)p"p^/p'^ and 
c{z;0 rdwxiw^) 



z-^ ' Jo WW 



& c(0;O = l (A.2) 



All non-universal information of the cut-off function is stored in the residues of the function c 
for 3fJz > 2. These residues are, however, of no relevance for universal properties of the system 
and further description of them is needed. The scope of this appendix is to apply the Mellin 
transform to the evaluation of the logarithmic asymptotics of the integrals introduced in the main 
text. The interested reader is referred to [12] for more details on this technique. 

Appendix A.7. Evaluation ofV^^l ^ 

After taking the Mellin transform with respect to the two cut-off scales m and x t (44) be- 
comes 



9(0) 

(rfe.);2,C 



Dom'pid-l)c{z;0) f dS 1 f dC r(-C) (>^r)« 



J (2 7r 



(A.3) 



xz J {2Try \qd-+^ J {2m) qd+^-^C 

The momentum integral can be performed for 

^{z-2C) < & ^z < (A.4) 

The asymptotics for 



XT 



« 1 (A.5) 



corresponds to a shift to the left of the contour of the inverse Mellin transform in In order to 
extricate the leading order of the asymptotics it is sufficient to approximate the inverse transform 
with the residue of the first pole which is encountered for — z/2 (3Rz < 0): 



^ ,1 (^.Y r r(|)r(i j) _ r(i-j) 



Finally the limit of vanishing m requires shifting the contour of the inverse MelUn transform in 
z towards positive values of 3? z. The result is 

V(°^°) {x,t) = -^hi—+ sub-leading (A.7) 

2 XT 
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(2) 

Appendix A.2. Evaluation of V^^^ ^ 
Proceeding as above (36) becomes 

i (2 7r)''' 1^ J i2m) q-^i ] g<^+^+2 "^^^"^-^ ^^'^^ 
The momentum integral exists for 

SRC = -2 + 2 < & < -2 (A.9) 

The asymptotic expression of the logarithmic part of the vertex in the limit (A.5) is obtained by 

approximating the inverse Mellin transform in C with the residues of the poles for SR^ = (2:+2)/2 
and 5RC = z/2. The resulting pole in z is double with residue 

^(?r.) ,Xj,t) = ->ct\a. Vrn?xt d^^ ■ d^^ 

+ %~/2n/-l) '^"'^^''''^' ^ "^^^-^^^^^ + sub-leading (A.IO) 

with 

d+ z(l — p) — 1 x-' 
Note that whenever acting on translational invariant functions 

^(tr.) ,Xj,t) = Ktln Vm?7tt An 

+ E'^"'('^^^'^'^)1^^4^-?^c.f +^"b-leading (A.12) 

Appendix B. Deiuiition of quadratic Casimir 

Let i be the particle label i = 1, . . . ,n and a, /3 be index of vector components of individual 
particles in R'^. Upon defining 

Qf^:=xfd^, & ti:=Qf^ (B.l) 

i 

the dilation operator acting on the n-particle space M.^ ^ is 

71 

i=l 

On the same space, the representation of the quadratic Casimir of SU{d) is 



2J »J 
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Using the dilation operator, this latter can be also written in terms of the quadratic Casimir of 

SU{n-l): 

4^2) = 4^(i-i) + l^^TTiy^" + din - 1)] (B.4) 
The generators of SO{d) for the i-th particle in d dimensions are 

17^ ■■=97^-97^ (B.5) 

The generators of SO{d) for n particles are additive functions of the generators of SO{d) for 
single particle: 

N 

K^=Y,lt^ (B.6) 

i=l 

The quadratic Casimir is 

1 " 

4oS) ■= 2^"^^-0 = E {t-^ (3i9j)} - tr (flifl*) (B.7) 
A excellent introduction to group theoretic methods for many-body systems can be found in [ 1 5 J . 



Appendix C. Derivation of (38) 

Omitting the Laplacian term as it vanishes on harmonic polynomials [V^^^^^^ ^] reduces to 



d + 2(l-p)-l 
2(d + 2) (d-l) 



(C.1) 



Unfolding the differences and dropping differential terms vanishing on trans- 

lation invariant functions, yields 



E 



2(l-rfp) 
rf + 2(l-p)- 



2(1 -dp) xfx^ + xfx^^ 



- ^ l^Y' ^ rf + 2(l-p)-l 

Applying the definitions of appendix Appendix B, finally gives 

((i+l~2p) 



•^-i cr. 



(C.2) 



{d+2) 



(2,„) _ d+_(d--2)_p (2,„) (d- l)(rf+2)p ^ _ 
^so{d) rf+i_2p^s^^W rf(rf+l-2p) """^ 

The equality (B.4) permits to express this result in terms of C^sl^f^^-i) ■ 



(C.3) 
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